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In integrable many-particle systems, it is widely believed that the stationary state reached at late 
times after a quantum quench can be described by a generalized Gibbs ensemble (GGE) constructed 
from their extensive number of conserved charges. A crucial issue is then to identify a complete set 
of these charges, enabling the GGE to provide exact steady state predictions. Here we solve this 
long-standing problem for the case of the spin-1/2 Heisenberg chain by explicitly constructing a GGE 
which uniquely Hxes the macrostate describing the stationary behaviour after a general quantum 
quench. A crucial ingredient in our method, which readily generalizes to other integrable models, 
are recently discovered quasi-local charges. As a test, we reproduce the exact post-quench steady 
state of the Neel quench problem obtained previously by means of the Quench Action method. 

PACS numbers: 02.30.Ik,05.70.Ln,75.10.Jm 


Introduction. Understanding and describing the equi¬ 
libration of isolated many-particle systems is one of the 
main current challenges of quantum physics. The pres¬ 
ence of higher conserved charges (above the Hamiltonian) 
is linked to the absence of full relaxation to a thermalized 
state; the conjectured appropriate framework to charac¬ 
terize the steady state properties in such a situation is 
the Generalized Gibbs Ensemble (GGE) [1], in which all 
available charges are ascribed an individual ‘chemical po¬ 
tential’ set by the initial conditions, and the steady state 
is the maximal entropy state fulfilling all the constraints 
associated to the conserved charges [mz]. The basic 
idea underlying the GGE is as follows. Let H ~ be 
the Hamiltonian of an integrable model, and a set 

of conserved charges fulfilling = 0. The sit¬ 

uation we are interested in is that of a quantum quench, 
where we initially prepare our system in the ground state 
1111(0)) of a local Hamiltonian Hq and then consider uni¬ 
tary time evolution with respect to our integrable Hamil¬ 
tonian 


|vl/(t)) = e-^‘|vk(0)). (1) 

We assume that we are dealing with a generic case, where 
in the thermodynamic limit |1E'(0)) cannot be expressed 
as a linear combination of any finite number of eigen¬ 
states of H. At late times after the quench expectation 
values of local operators approach stationary values 

(0)^= lim(vE-(t)|0|vI/(t)). (2) 

t—^oo 

The GGE hypothesis asserts that these expectation val¬ 
ues can be calculated as {O)^ = Tr(gGGEO) from a sta¬ 
tistical ensemble with a density matrix 


PGGE = ^ exp 




(3) 


Here Z is a. normalization, and the Lagrange multipliers 
fin are fixed by the initial conditions 


Tr(eGGELf(”^) ,. («'(0)|7J(")|^(0)) 

limth --VT-- = hmth — 


N 


N 


(4) 


where N is the system size and limth denotes the ther¬ 
modynamic limit N —> oo. Eqns. ([^ are a direct con¬ 
sequence of the fact that are conserved charges. 

While the GGE hypothesis has been successfully verified 
for many systems mappable to free particles, in interact¬ 
ing theories such as the spin-1/2 Heisenberg XXZ chain 
the question arises, precisely which charges need to he in¬ 
cluded in In Refs. |17L 1181 120] a GGE based on 
the known conserved local charges |25] was constructed 
and used to determine steady-state averages of observ¬ 
ables [50] ■ Subsequent analyses Emm by the Quench 
Action (QA) approach [T3] demonstrated that this GGE 
fails to predict the correct steady state properties. This 
failure was shown to be related to the existence of bound 
states [U] (see also ESJ [30] ), which are known to be a 
generic feature in quantum integrable models. These re¬ 
sults posed the question whether the GGE is conceptu¬ 
ally faulty, or whether there could exist hitherto unknown 
charges that need to be taken into account in its construc¬ 
tion. 

In this Letter, we settle this issue by explicitly showing 
how to repair the GGE in Heisenberg chains, by com¬ 
plementing it with recently-discovered additional fami¬ 
lies of conserved charges [3T]- Grucially, these “quasi¬ 
local” charges fulfil a weaker form of locality than the 
previously known ones. We derive a set of fundamen¬ 
tal identities between the initial-state expectation values 
of these charges, and the density functions characteriz¬ 
ing the steady state. An explicit test of our construction 
is provided by a quantum quench from the Neel state 
to the XXZ chain: we demonstrate that our GGE cor- 
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rectly recovers the stationary state, the form of which is 
known exactly from the QA approach [ 2 TJ [32] . In this 
way we completely resolve the above-mentioned conun¬ 
drum. Our construction shows that quasi-local conserved 
charges are in fact crucial for understanding the non¬ 
equilibrium dynamics of quantum integrable models. 

Anisotropic spin-1/2 Heisenberg chain. We shall con¬ 
sider a completely generic quench protocol from an initial 
pure wavefunction I'I'o), which is unitarily time evolved 
according to the Hamiltonian 

^ = - 1 )] ■ (5) 

i=i 


Here J > 0, cr“, a = cc, y, z are Pauli matrices acting on 
spin- 1/2 degrees of freedom, and we consider anisotropy 
values in the regime A = cosh ( 77 ) > 1. The Hamiltonian 
([^ can be diagonalized by Bethe Ansatz [331131] • Im¬ 
posing periodic boundary conditions, energy eigenstates 
I A) with magnetization H — M are labeled by 

a set of rapidities A = {Xk}^:^i satisfying the Bethe 

sin(Aj—Afc-l-ir;) . _ 

sin(\j-X^-iri) ’ J ~ 


equations = - 


= -n“i 


1,..., M. The momentum and energy of a Bethe state 
are Px = = E^ie(Aj) where p{X) = 


i In 


sin(A—277/2) 
sin(A+ 277 / 2 ) 


and e(A) = — j 7 rsinh( 77 )ai(A), where 


^n(A) — 


2 sinh(7W7) 


27r cosh(n 77 ) — cos (2A) ’ 


( 6 ) 


the form (§, Q with charges Q. The initial val¬ 
ues = limth IV~^(^'(0)|i7*'”^|'I'(0)) of the conserved 
charges are conveniently encoded in the generating func¬ 
tion [T3] 

= limth f) dxT{\ + f) Ivko), 

\ k 

fc =0 

Given the GGE density matrix, a “microcanonical” de¬ 
scription of the steady state can be obtained by perform¬ 
ing a generalized Thermodynamic Bethe Ansatz (GTBA) 
[111137], see SM for a brief summary. This results in a 
representative eigenstate Ip^ge) labeled by root density 
functions Pqqe’ which has the property that for any local 
operator O 


Tr(C>pGGE) — (PggeI^Ipgge)- (9) 

Within the GTBA formalism macrostates can be de¬ 
scribed either by root densities of particles, or by den¬ 
sities of holes. Holes can be, loosely speaking, under¬ 
stood as analogues of unoccupied states in models of 
free fermions. In terms of the latter the state IPqge) 
is parametrized in terms of the set of positive functions 
{p!.u. In mm it was found that the initial data (j^ 
directly determines the hole density of I-strings (i.e. va¬ 
cancies of unbound states), according to the remarkable 
identity 


Solutions A to the Bethe equations are closed under com¬ 
plex conjugation and consist of so-called strings A)J’“ = 
AS + f (" + 1 - 2 a) + a = 1 ,... ,n and AS e M. 

Here index a enumerates a string, n is the string length, 
a counts rapidities inside a given string and deviations 
i5S’“ are (for the majority of states) exponentially small 
in system size [TS] [33] [33] . The string centers AS lie in 
the interval [— 7 r/ 2 , 7 r/ 2 ). In the thermodynamic limit 
IV —>■ c» with M/N fixed one can describe a state not 
in terms of individual rapidities, but rather in terms of a 
set of functions p = {pn}//! representing string densities 
(see Supplementary Material (SM) for more info). 

Ultra-local GGE treatment. Exactly-solvable Hamil¬ 
tonians such as ([^ can be embedded [35] in a commuting 
family [T(A),T(A')] = 0 of transfer matrices (defined in 
@). The Hamiltonian and an infinite number of mutu¬ 
ally commuting conserved charges are obtained by 

=-d-xlnT{-tX) (7) 

n\ A=^ 

with the Hamiltonian reading H = iUlMllljjU) ^ These 
charges are ultra-local in the sense that they can be writ¬ 
ten as Tfl™) = where the operators 

act nontrivially on a block of at most m sites adja¬ 
cent to j. The GGE constructed in m IB] was of 


plliX) = a/X) + ^ [H-^nA + f) + H'^«(A - f)] . 

( 10 ) 

All other hole densities are fixed by the maximum en¬ 
tropy principle under the constraints Q. 

Quench Action treatment. The above GGE treatment 
should be compared to an independent calculation us¬ 
ing the QA method |I4j . For a generic quench problem, 
given an initial state l^/o), the time-dependent expec¬ 
tation value of a generic local observable O can be ex¬ 
pressed as a double Hilbert space summation 

(4'(t)|C>|4'(t)) = ^ e-^^-'^^'e*(“^-“^')‘(A|C>|A'), 

A.A' 

( 11 ) 

where Sx = — ln(A|tI>o). Here, |A) are eigenstates of 
the Hamiltonian driving the post-quench time evolution. 
Exploiting the fact that in the themodynamic limit, the 
summation over eigenstates can be written as a func¬ 
tional integral over root densities, which can be evalu¬ 
ated in a saddle-point approximation (becoming exact in 
the thermodynamic limit), one finds in particular that 
the steady-state expectation values of observables a long 
time after the quench can be obtained as 

hm limth i'fiQl O |vk(t)) = {pIi\0\pIi). ( 12 ) 
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Here |pqa) is an eigenstate minimizing the QA S'qaIp] = 
2S[p] — 5yy[p], where S'[p] = limth Re S\ is the extensive 
real part of the overlap coefficient in the thermodynamic 
limit and 5'yy[p] is the Yang-Yang entropy of the state 
PSI1551 [55] , For the Neel to XXZ quench, the exact 
overlaps were obtained in [35] and used in EH [3H to 
obtain the exact saddle-point densities represent¬ 
ing the steady state. Crucially, one finds [211 Eg that 
Pgge ^ Pqa^ which in turn leads to different predic¬ 
tions for physical properties such as spin-spin correla¬ 
tors. This demonstrated that the ultra-local GGE does 
not correctly describe the steady state after a generic 
quantum quench in the XXZ chain. 

Constructing a “quasi-local” GGE. Very recently |3T] 
(see also |53H33] 1 hitherto unknown conserved charges 
of the isotropic (A = 1) Heisenberg model were discov¬ 
ered. These operators are not local in the sense that they 
cannot be represented as a spatially homogeneous sum 
of finitely supported densities, but rather quasi-local, 
meaning [3T] that their Hilbert-Schmidt norms scale lin¬ 
early with system size and their overlaps with locally- 
supported operators become independent of N in the 
limit of large system size. Moreover, they are linearly in¬ 
dependent from the known local charges generated from 
the spin-1/2 transfer matrix. Until now, the impact of 
these charges on local physical observables has not been 
quantified. 

Our first step is to construct a family of quasi-local con¬ 
served charges for A > 1 by generalizing the procedure 
of [31]. The starting point is the q-deformed L-operator, 


L{z, s) = — ( sinh(z) cosh {rjs^) 0 (13) 

smh r] V 

-l-cosh(z) sinh (rjs^) 0 smh(-g)(s~ 0 a^-\- s"*" 0 , 


whose auxiliary-space components are given by q- 
deformed spin-s representations with s = |,1,|,..., 
obeying commutation relations [s+,s“] = [ 2 s^]g, 

= ±s^ and acting in a ( 2 s -|- l)-dimensional 
irreducible representation Vs = = lsp{|fc);A: = 

-s, ...,s}, 

s'" \k) = k\k), 3=*= |fc) = ^[s -f 1 ± A:]q[s =F k]q |fc ± 1), 

(14) 

with [x\q = sinh( 77 x)/sinh( 77 ). By means of higher-spin 
auxiliary (fused) transfer matrices defined via ordered 
products of L-operators 


Ts{z) = Tra [La,i{z, s)... La,N[z, s)], (15) 

(where Ti/ 2 [z) = T{z) was used in Q) we define families 
of conserved operators 

As(A) = T-\X){Uz^)niz+)l 4 = ±1 + (16) 

with T'g{z) = dzTs{z) and {•} denoting the traceless part. 
The normalization reads rs(A) = f{—is-\-^)r]-\-iX)f{{s-\- 


^)r]-\-iX) with f(z) = (sinh(z)/sinh( 77 ))^. In ^ it was 
shown for the isotropic case that these charges are quasi¬ 
local for all s = 5 ,1, |,... and A S M. A rigorous proof 
for general A > 1 is currently under construction [35] . 

A central piece of our work is the extraction of the ther¬ 
modynamically leading part of the quasi-local charges 
{As}^i /2 when operating on an arbitrary Bethe state. 
It proves useful to resort to the so-called fusion relations 
[35H33] (T-system) for higher-spin transfer matrices, 

Uz + f )r,(z - f) = f{z + {s+ ^)r,)f{z -{s+ i), 7 ) 

+ Te-i/2{z)Ts.^.l/2{z), (17) 

with the initial condition Tq(z) = f{z). There exists a 
closed-form solution to the above recurrence relation in 
terms of Baxter’s Q-operators [47] 

Ts{z) = Q{z -\-{s-\- ^)v)Qiz - (s -f 4)??) (18) 

X /(z + (€ — s)ri) 

^ Q{z + (e - s - f)r]jQiz -^ie-s + 4)V)' 


The eigenvalues of the Q-operators (in what follows, in 
view of commutations [Ts{zi),Q{z 2 )] = 0 Vs,Zi € C, we 
slightly abuse notation by using the same symbol for an 
operator and its eigenvalue) are determined by the posi¬ 
tion of Bethe roots Q{z) = Jlfeii (z + *Afc). A key 
observation is that, in the thermodynamic limit, the spin- 
s transfer matrix evaluated at z^ (z/") is simply given by 
the i = 0 {£ = 2s) term in the sum in Eq. (18). This 
then gives 


limth Ts(z)l') = limth/(±(s+ 5)77+ *A) 


QiTsr] iX) 


Q(±sr7 - 


- iX)' 
(19) 


The latter analysis is consistent with 
limth'rU^(A)Ts(z/')Ts(z/") = 1 , representing a ther¬ 
modynamic version of an inversion identity (see |31j l 
that can be proven in an entirely operatorial way, 
without making reference to the Bethe eigenstates. At 
this point it is convenient to define modified conserved 
operators 


A«(A):=ri-)(z)/)Ti+)'(z+), (20) 


where (z) have the same structure as (15) but involve 
L-operators L(=*=)(z, s) = L(z, s) sinh ( 77 )/[sinh (z ± 577 )]. 
In thermodynamic limit a quasi-local conserved op¬ 
erator As (A) only differs from As(A) 
of identity, As(A) = As(A) -I- ts(A)]l, 

2 s sinh(( 2 a-H)i 7 ) -l/'v', ^ 

2 S-I -1 sinh^ (77) 

parameter family of conserved charges 


by a multiple 
with ts(A) = 
define a two- 


1 


F("+i) = -d^XsiX) 


A=0 


( 21 ) 


By construction we have [Hi^\ = 0 and 

precisely recover the ultra-local conservation laws ([^. 
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We are thus in a position to define the density matrix 
of our GGE. It is given by 


1 

9GGE = exp 


- 5] 


/2 


( 22 ) 


where the Lagrange multipliers /3® are fixed by initial 
conditions of the form Our assertion is that ( |2^ 
provides a correct description of the stationary behaviour 
after a general quench to the spin-1/2 XXZ chain (in 
the regime A > 1). In order to prove this it suffices to 
establish that the initial values of our conserved charges 
uniquely specify a macrostate. 

Let us now derive the main result of our Letter. Anal¬ 
ogously to what was found in for the ultra-local 

charges, the values of the quasi-local charges associated 
with a spin-s transfer matrix are in 1-to-l correspondence 
with functions P 2 shW^ '"^hich in turn specify (see SM) a 
unique macrostate (namely the GGE saddle-point state). 

Our starting point is the following expression for the 
spectrum of {As}/li/ 2 ) valid for large system size (cf. 
Eq. 


Xs{X) = -idxlog 


M 


= E 

k=l 


Q{-st] -I- iX) ^ 
Qisij + iX) 

2 sinh ( 2 s? 7 ) 


o{N) 


(23) 


cos(2(Afc -I- A)) — cosh(2sr?) 


+ o{N), 


Starting from Eq. (23), working in the thermodynamic 


limit under the string hypothesis and making use of Bethe 
equations, one arrives at (see SM) 

= «2«(A) + ^ [Ofo(A + f) + ^I^(A - f)] , 

(24) 


where s = 5,1 ,|,— The right-hand side of (24) is 
determined by the expectation values of the quasniocal 
charges on the initial state. 


n'f°{X) = limth 


(4'o|X,(A)|4'o) 

N 


(25) 


This is a generalization of Eq. (10) to arbitrary spin. 


Note the remarkable fact that this correspondence is valid 
for a generic initial state l^/o). As a consequence, the 
family of quasi-local charges {^s}“ 1/2 completely deter¬ 
mines the postquench stationary state through the GGE 
and gives the latter’s predictions identical to those com¬ 
ing from the QA. 

Neel quench. An explicit example of our construction 
if provided by the quench from the Neel state 

i4/o) = ^ (inn. •.) + inn •■•))• ( 26 ) 

Here the root distributions characterizing the stationary 
state have been previously determined by a QA calcula¬ 
tion HUES]. In order to demonstrate that our GGE re¬ 
covers these known results we need to compute the gener¬ 


ating functions (25). Here we can repeat the logical steps 
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FIG. 1. Comparison of methods: QA method versus im¬ 
proved GGE predictions. Colored lines pertain to the re¬ 
fined GGE calculation with systematic addition of higher- 
spin families of quasi-local charges for the local cor¬ 

relation function (crfcrl) in the regime A > 1 (left panel). 
Labels in GGEs indicate the maximal auxiliary spin s for 
the charges being included in the GGE computation. 

The right panel displays the relative differences S(a-}rT§) = 
({cio- 3 )QA-(crfcr|)GGBg)/{fric- 3 >Q 2 i in logarithmic scale. We 
used the mapping between correlation functions and the set 
of densities p given in |50| . 


from the calculation for s = 1/2 in [Tallin] by studying 
the spectrum of associated auxiliary transfer matrices. 
This calculation can be found in the SM. Substituting 
the results obtained in this way into (24) gives perfect 
agreement with the known QA results. 

Towards a truncated GGE. In m it was argued that 
for the purpose of describing finite subsystems in the 
thermodynamic limit ultra-local GGEs can be truncated 
by retaining only a finite number of the “most local” con¬ 
served charges. An obvious question is whether a similar 
logic can be applied to our quasi-local GGE. As a first 
step towards understanding this issue, we have calculated 
the next-nearest spin correlation function in the steady 
state after a Neel-to-XXZ quench for several GGEs trun¬ 
cated in the s direction. In Fig. [^we show the results of 
these calculations for A > I. The data clearly shows that 
adding subsequent families of quasi-local charges results 
in a rapid convergence of the corresponding truncated 
GGE result to the exact value. 


Gonclusions. We have shown how to construct an ex¬ 
act GGE describing the stationary state after generic 
quantum quenches to the spin-1/2 Heisenberg XXZ 
chain. Our GGE is built from an extended set of lo¬ 
cal and quasi-local charges. We have shown that our 
construction resolves previously observed discrepancies 
between predictions for steady state expectation values 
by an exact QA treatment on the one hand, and a GGE 
restricted to ultra-local charges obtained from the trans¬ 
fer matrix of the spin-1/2 chain on the other hand. Our 
results provide unambiguous proof that the recently dis¬ 
covered quasi-local charges have a non-negligible impact 
on the relaxation processes of strongly-interacting many- 
body quantum systems in one dimension. 
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We begin by recalling some of the fundamental equa¬ 
tions of the Bethe Ansatz solution of the XXZ chain, 
together with few results obtained in EHEa which are 
used in the bulk of our paper. 


Thermodynamic limit of Bethe equations 

In the thermodynamic limit, the Bethe equations for 
the XX Z chain read [33 [351 [31] 


^ ^ {^nm * t (^t) 

m—1 

for n > 1, where p„,t(A) = /o„(A) -I- Pn,hW and Pn,Pn,h 
are respectively the particle and hole densities of n- 
strings. The convolution is defined by (/ * g) (A) = 
I-i% ~ dip)- The kernels are 

^nm(^) ~ (1 ^nm)^\n—m\i^) “f ‘^^\n—m\+2i^) T ■ ■ ■ 

-j- 2Q.^_|_m_2 (A) -f Uji +miX), (S2) 

with an defined in in the main part of the manuscript. 
A convenient rewriting is in the decoupled form |36j 

Pn(l + ?7«) = S * iVn-lPn-l + Vn+lPn+l), (S3) 


GTBA for the GGE 


The generalized TBA for the GGE based on local 
charges proceeds as a standard TBA, but now with the 
effect of additional charges beyond the Hamiltonian being 
taken into account by additional parameters /3„ (chemi¬ 
cal potentials) in the GGE density matrix. By applying 
the standard maximal entropy reasoning using these con¬ 
straints results in the GTBA equations Eiiisa 

ln(? 7 „) = -(5„,i(s*d)-|-s*[ln(l-|-?7„_i)-|-ln(l-|-?7„+i)] , 

(S6) 

for n > 1, where r]o{X) = 0 and s(A) is defined in ( j^ . 
The driving term originating from ultra-local charges is 
remarkably only present in the first integral equation and 
is specified by the chemical potentials Pm, m> 2, 

OO 

d{X) = ^ ^ Pm . (S7) 

m—2 


As shown in miisi], the ultra-local charges associated 
to the spin-1/2 transfer matrix completely fix the density 
of holes of 1-strings to but leave all higher hole 
density functions Pnh, with n > 2, undetermined. As 
explained in detail in [33] , the GTBA system of (S6) and 
(S3) for the GGE can then be solved (using the constraint 
Pi h = pf°fi to eliminate the unknown driving term in 

m). 


for n > 1, where rjn = Pn,h/Pn- The A-dependence is 
left implicit and we use the conventions rjoiX) = 1 and 
po(A) = (5(A). The kernel in ( [S^ ) reads 


s(A) 





g—2i/cA 
cosh(fc77) 


(S4) 


The set p = {Pn}))Ti represents an ensemble of states 
with Yang-Yang entropy 


GTBA for the Quench Action 

In the case of the QA treatment, the GTBA equations 
take the form 

ln(p„) =dn+S* [ln(l -h T]n-l) + ln(l -f Pn+l)] , (S8) 

where n > 1. The driving terms are given by the exact 
overlaps of Bethe states with the initial state. In the 
specific case of the Neel quench, these are given by 


OO pTr/2 

Syy [p]=Ny 2 dA [pn,t (A) In p„,* (A) 

/2 

-pniX)llipniX) - Pn,hiX)lnpn,hiX)] . (S5) 


dniX) = Y. 


-2*fcA tanh(fc?7) 


[(-ir-(-i)l 


feGZ 


= (-!)" In 


rdI(A)1 


r^i(A)i 

mx)\ 


mx)\ 


(S9) 


An important point to bear in mind is that the Bethe 
equations (S31 relate the set of densities p to the set of 
hole densities pf^. Knowing one of these two sets is thus 
sufficient to completely determine a given state. This 
point is crucial to understand the effects of constraints 
coming from ultra-local and quasi-local charges, as is ex¬ 
plained below. 


where -dj, j = I,...,4, are Jacobi’s "d-functions with 
nome 

Note the difference between this GTBA and the one as¬ 
sociated to the ultra-local GGE: from the exact QA treat¬ 
ment one obtains GTBA equations with driving terms at 
all string lengths n, as a result yielding a different set of 
the steady state densities. 
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Relating quasi-local charges and hole densities 


terms p2s,h- 


In the main text in (23) we provided the expectation 


values of the quasi-local charges on Bethe eigenstates in 
the limit of large system size. Strictly in the thermody¬ 
namic limit one obtains 


limth ^ (A|Xs(^) I A) 

. M 




2 sinh ( 2 s? 7 ) 


N cosh (2(zfe — ip)) — cosh (2s?7) 

min(n,2s) 


(SIO) 


oo .7r/2 mm(n, 2 s) 

= -27r^/ dXpniX) ^ a|„_2s|_i+2j(A-f/x), 

n^lJ-TT/2 

(Sll) 

with an defined in (§. In the last equality we accounted 
for Zk = —iXk, used the string hypothesis and the fact 
that the expectation values can be written as 

(A|i,(M)|A) 

sm(Afc + p- isp) J 

Using conventions for the Fourier transform 

fK/2 


m = l_ 


-7r/2 


dAe^*'=-'/(A) , kGZ, 




— 2ik\ 


m, Ae[-f,f) 


feGZ 


one can map Eq. (Sll) to Fourier space 


limth ^ {\\Xsip) I A) 


(513) 

(514) 

(515) 


OO min(n,2s) 

= ^-\k\r,(\n-2s\-l+2j) ^ 

/cGZ n—1 j — 1 


using that a^ik) = By performing the sum over 

j and using that \n — 2s| -I- 2min(n, 2s) = n -|- 2s, one 
finds 


limth ^ WXs{p) |A) 


(S16) 


= E 


g—z2fc/j. 


feGZ sinh(|fc| 77 ) 


OO 

(' 


^-\k\rj{n+2s) _ g-|/c|r;|n-2s| 


Using the thermodynamic Bethe equations (cf. (S3)) in 
Fourier space, 


PnAk) = 


1 


2 cosh(fc77) 


(/3„-i.h(fc) +/On-hi.h(fc)), (S17) 


where po^h{k) = 1, one can observe a cancellation of all 
terms with an exception of an expression given solely in 


limth ^ A\Xs{p) I A) 


= E 


^—i2k^ 

cosh.{kri) 


{p2s,h{k) 


_ p-2s|fe|r; 


)■ 


(S18) 


The quasi-local conservation laws make the left-hand side 


of Eq. (S18) equal to the generating function of the 


charges on the initial state as stated previously in (25), 
leading to 

Uf“(fc)C0®h(M=p2..h(A:)-e-2«l'=IU (S19) 

TT 

Taking the inverse Eourier transform produces the main 
result of our Letter, namely the identification given by 


(24). 


Truncated GGE 

Eor practical reasons it is useful to determine a GGE 
ensemble by including only a finite number s of quasi¬ 


local charges Using Eq. (24) for s = l,...,s 

we fix the distributions of holes pn,h = p'Ah string of 
lengths n = 1,..., 2s. These restrictions can be in turn 
used as a driving term for the following GTBA equations 
(analogously to what has been done for the case s = 1/2 
in [1111231^) 

ln(77„) = s * [ln(l -I- 77„_i) -b ln(l -b 77„+i)] n > 2s , 

(S20) 

with 


(S21) 


'Po 

Vn = — n<2s, 

Pn 


where the functions pn solve the Bethe equations 
Eor any s > 1/2 we then apply an iterative procedure to 
find a solution for all the all the 

ultimately leading to the results shown in figure 


Neel initial state 

As the Neel state is a simple product state we can 
evaluate all scalar products in the “quantum” spaces 
pertaining to the physical spin-1/2 degrees of free¬ 
dom. This leaves us with a staggered product of diag¬ 
onal components of an auxiliary two-channel L-matrix 
Ls ( 21 ,^ 2 ) = (z 2 ,s), producing a transfer 

matrix Ts(zi,Z 2 ) operating on two copies of auxiliary 
spin-s spaces Vs G Vs, 


^s{zi,Z2) = L/^(2;i,Z2)L/'^(zi,Z2). 


(S22) 


This is a local unit of the ‘boundary partition function’ 
Zs{zi,Z 2 ) = limth lV“^TraTs(zi, 2 : 2 )^^^, whose large-A^ 















limit on I? := {( 2 :^, 2 ;^); A S K} C is dominated 
by a non-degenerate unit eigenvalue Ag {z^ ,z^) = 1 of 
Ts(zi, 2 ; 2 )|-d, implying 

= 5,,Zg(zi,z2)b = ^KAg(zi, 22 )]p. (S23) 

It also helps noticing that Tg(zi,Z 2 ) enjoys a C/(l)- 
symmetry, with the leading eigenvalue always residing in 
the largest (2s -I- 1)-dimensional subspace. Closed-form 
results can be readily obtained in the cases s = ^, 1 


suppressing for higher spins s > 3/2 which be¬ 

come quickly cumbersome expressions. For practical im¬ 
plementation of the truncated GGE it however suffices 
to evaluate them numerically. For a class of initial states 
which are given in the Matrix Product State form this 
can be done efficiently by e.g. employing the method 
outlined in pT)] . 




— sinh ( 277 ) 

1 — 2 cos (2A) -I- cosh {2r]) ’ 

2 sinh ( 3 ? 7 ) 

3 cos (2A) — cosh ( 77 ) — 2 cosh ( 3 ? 7 ) 


(524) 

(525) 




